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Abstract 

Using pointwise semigroup techniques of Zumbrun-Howard and 
Mascia-Zumbrun, we obtain sharp global pointwise Green function 
bounds for noncharacteristic boundary layers of arbitrary amplitude. 
These estimates allow us to analyze linearized and nonlinearized sta- 
bility of noncharacteristic boundary layers of one-dimensional systems 
of conservation laws, showing that both are equivalent to a numerically 
checkable Evans function condition. Our results extend to the large- 
amplitude case results obtained for small amplitudes by Matsumura, 
Nishihara and others using energy estimates. 

1 Introduction 

Boundary layers appear in many physical settings, such as gas dynamics, 
MHD, and rotating fluids; see, for example, the physical discussion in |SGKOj . 
In this paper, we study the stability of boundary layers assuming that the 
boundary layer solution is noncharacteristic which means that signals are 
transmitted into or out of but not along the boundary. Specifically, we on- 
sider a boundary layer, or stationary solution, 

(1.1) u = u{x), lim u{z) = u+, u{0) = uo 
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of a system of conservation laws on the quarter-plane 



(1.2) 



ut + = {B{u)u^) 



xA > 



u, f & R", B G M"^", with initial data u{x, 0) = g{x) and Dirichlet boundary 
condition 



A fundamental question is whether or not such boundary layer solutions are 
stable in the sense of PDE, i.e., whether or not a sufficiently small pertur- 
bation of u remains close to u, or converges time-asymptotically to u, under 
the evolution of (11. 2p . 

Long-time stability of boundary layers has been considered for scalar 
equations in [LN[ ILYj and for the equations of isentropic gas dynamics in 
pNl IKNZ] . The latter results, obtained by energy estimates, apply to ar- 
bitrary amplitude layers of "expansive inflow" type analogous to rarefaction 
waves, but only to small-amplitude layers of "compressive inflow or outflow" 
layers analogous to shock waves or "expansive outflow" type. For general 
symmetric hyperbolic-parabolic systems, stability of small-amplitude non- 
characteristic boundary layers has been shown in multi-dimensions for strictly 
parabolic systems in |GGj . and in one dimension for partially parabolic ("real 
viscosity") systems in [Rlj . Here, in the spirit of results obtained for shock 
waves in |ZH[ IMZ3[ IMZ4j , we show for general strictly parabolic systems of 
conservation laws that linearized and nonlinear stability are equivalent to a 
generalized spectral stability condition phrased in terms of the Evans func- 
tion associated with the linearized equations about the wave, independent of 
the amplitude of the boundary layer in question. 

The Evans condition is readily checkable numerically, and in some cases 
analytically; see [Eni [Br2l EEZl [BDUl MuZi IBHRZl [HLZ] . In particular, 
stability of small-amplitude uniformly noncharacteristic boundary layers has 
been shown for general hyperbolic-parabolic systems in multi-dimensions in 
[GMWZl] using elementary Evans function arguments (convergence to the 
constant layer). An exhaustive numerical study for isentropic gas layers 
in one dimension has been carried out in |CHNZj . with the conclusion of 
stability for arbitrary amplitudes. 

Our method of analysis is by pointwise Green function methods like those 
used in |ZHt IMZ3t IMZ4] . and especially |HZ] . to analyze the stability of 



(1.3) 



u{0,t) = h(t). 
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viscous shock layers. Similar results have been obtained for the related small- 
viscosity-limit problem in [GRl IMeZl IGMWZ2j . In particular, we point to 
the analysis of Grenier and Rousset |GRj as using pointwise Green function 
estimates very similar to those that we use here, though adapted for different 
purposes. 

1.1 Equations and assumptions 

Consider a viscous boundary layer, a standing-wave solution (11. ip of a gen- 
eral parabolic system of conservation laws fll.2p . Assume, similarly as in the 
treatment of the viscous shock case in |HZj : 

(HO) f,BeC^. 
(HI) Rea{B)>0. 

(H2) cr(/'(M+)) real, distinct, and nonzero. 

(H3) Rea{-ikf{u+) - k'^B{u+)) < -9k^ for all real k, some ^ > 0. 
(H4) Solution u is unique. 
Here, (H2)(iii) corresponds to noncharacteristicity. 



1.2 Linearized stability and the Evans function 

After linearizing II. 21 about the stationary solution u, we obtain the linearized 
equation 

(1.4) ut = Lu := -{Au)^ + (Bu^)^, A,B e C^, 
where 

(1.5) B := B{u) 
and 

(1.6) Au := dF{u)u- dB{u){u,u^). 

Definition 1.1. The boundary layer u is said to be linearly asymptotically 
stable, if u{-,t) approaches as t ^ oo, for any solution u of ( [i.^[ ) with 
initial data bounded in in some specified norm. 
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We define the following stability criterion, where D{X) described below, 
denotes the Evans function associated with the linearized operator L about 
the layer, an analytic function analogous to the characteristic polynomial of 
a finite-dimensional operator, whose zeroes away from the essential spectrum 
agree in location and multiplicity with the eigenvalues of L: 

(1.7) There exist no zeroes of D{-) in the nonstable half-plane ReX > 0. 

As discussed, e.g., in |R2j . under assumptions (H0)-(H4), this is equivalent 
to strong spectral stability, (j{L) C {ReX < 0}, (ii) transversality of u as a 
solution of the connection problem in the associated standing-wave ODE, 
and hyperbolic stability of an associated boundary value problem obtained by 
formal matched asymptotics. Here and elsewhere a denotes spectrum of a 
linearized operator or matrix. 

Our first main result is as follows. 

Theorem 1.2. Assuming (H0)-(H4), linearized asymptotic Ci —>■ 
stability, p > 1, is equivalent to fll.7p . 

Theorem 11.21 is obtained as a consequence of the following detailed, point- 
wise bounds on the Green function G{x, t; y) of the linearized evolution equa- 
tions (11.41) with homogeneous boundary conditions (more properly speaking, 
a distribution), defined by: 

(i) {dt — Lx)G = in the distributional sense, for all x,y,t > 0; 

(ii) G{x, t; y) — > 5{x — y) as t — > 0; 

(iii) ^(0, t- y) = 0, for all y,t>0. 
Denote by 

(1.8) a;/" < < ■ ■ ■ < a+ 

the eigenvalues of of the limiting convection matrix := df{u^). 

Then, our second main result is as follows. 
Theorem 1.3. Assuming (H0)-(H4) and stability condition (ll.7p . 



d2d^G{x,t;y) \ < Ce-''(l^-^l+*) 



n 



(1.9) 



+ C'(r 1-1/2 + |a|e-''l^l + |7|e-^l^l)(5;^ri/2g-0 



x-y-af. t)'^ /Mt 



k=l 



+ 




a+<0,a+>0 
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< |«|, I7I < 1, for some rj, C, M > 0, where denotes the posi- 
tive/negative part of X, indicator function X{\a't\>\y\} ^■^ Z^'" l^kA — Ivl 
and otherwise. 

1.3 Nonlinear stability 

Definition 1.4. The boundary layer u is said to be nonlinearly asymptotically 
stable if u{-,t) exists for all t > and approaches u as t —>■ 00, for any 
solution u of U.^) with initial data sufficiently close in some norm to the 
original layer u. 

Denoting by 

(1.10) < < ■ ■ ■ < 

the eigenvalues of of the limiting convection matrix := df{u+), define 

(1.11) e{x,t):= ^(l + t)-i/2e-k-4t|Vit^ 

a+>0 
a+>0 

and 

(1.13) ij^ix, t) := (1 - x{x, t))(l + |x - a+t| + t'^Y'^^ 

where x{x,t) = 1 for x G [0,a^t] and x{x,t) = otherwise and L > is a 
sufficiently large constant. For simplicity, take B identically constant. Then, 
our third and final main result is as follows. 

Theorem 1.5. Assuming (H0)-(H4), B = constant, and the linear stability 
condition fll.7p . the profile u is nonlinearly asymptotically stable with respect 
to perturbations g, h in initial and boundary data satisfying 

\g{x)\ < Eo{l + |x|)-3/^ mi < Eo{l + \t\r''\ \h'{t)\ < Eo{l + |t|)-^ 

for Eq sufficiently small. More precisely, 

(1.14) \u{x,t) -u{x)\ < CEo{9 + i;^+i;2)ix,t), 

where u denotes the solution of (11. 2p with initial data g = u + g and boundary 
data h = uq + h. 
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Remark 1.6. Pointwise bound fll.l4p yields as a corollary the sharp decay 
rate 



1.4 Discussion and open problems 

The case of boundary layers is quite analogous to the undercompressive shock 
case; in particular, pointwise estimates as in jHZj appear to be necessary 
to close the one-dimensional analysis by a linearized semigroup approach 
suitable for large-amplitude layers. (On the other hand, small-amplitude 
stability has been estabhshed using energy estimates in, e.g., |MNl IGG] .) A 
new feature of the present analysis as compared to those of |HZt IHRt IHRZj 
is the admission of perturbations in boundary as well as initial data. Open 
problems are extensions to systems with physical (partial) or quasilinear 
viscosity and to mult i- dimensional boundary layers. 

2 The Evans Function 

Before starting the analysis, we review the basic Evans function methods and 
gap/conjugation lemma. 

2.1 The gap /conjugation lemma 

Consider a family of first order ODE systems on the half-line: 



These systems of ODEs should be considered as a generalized eigenvalue 
equation, with A representing frequency. We assume that the boundary ma- 
trix B is analytic in A and that the coefficient matrix A is analytic in A as 
a function from Q into L°°(x), in x, and approaches exponentially to a 
limit A+(A) as x — > oo, with uniform exponentially decay estimates 



(1.15) 



u{x,t) - u{x)\lv < CEo{l + t)-^^'-p, 



1 < p < oo. 



(2.1) 



W = A{x, X)W, Xen and a; > 0, 
M{\)W = 0, Xen and x = 0. 



(2.2) 



\{d/dx)\A- A+)\ < Cie 



e\x\/C2 



for X > 0, < k < K, 
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Cj, 9 > 0, on compact subsets of f2. Now we can state a refinement of the 
"Gap Lemma" of [GZllKSj . relating solutions of the variable-coefficient ODE 
to the solutions of its constant-coefficient limiting equations 



(2.3) 



Z' = A+{\)Z 



as X — > +00. 

Lemma 2.1 (Conjugation Lemma [MeZj ). Under assumption f l2.2p . there 
exists locally to any given \q E Q a linear transformation P+(x, A) = I + 
0+(x, A) on X > 0, analytic in A as functions from Q to L°°[0, +oo), 
such that: 

(i) \P+\ and their inverses are uniformly hounded, with 



(2.4) \{d/d\y{d/dxfQ+\ < C{j)CiC2e'^\''\/^^ for x > 0, < k < K + 1, 



j > 0, where < 6 < 1 is an arbitrary fixed parameter, and C > and the 
size of the neighborhood of definition depend only on 9, j, the modulus of the 
entries of A at Xq, and the modulus of continuity of A on some neighborhood 
o/Ao e n. 

(a) The change of coordinates W := P+Z reduces fl2.ip on x > to the 
asymptotic constant- coefficient equations (12. 3p . Equivalently, solutions of 
(12. ip may be conveniently factorized as 



where Z^ are solutions of the constant- coefficient equations, and 9+ satisfy 
bounds. 

Proof. As described in |MZ3] , for j = k = this is a straightforward corollary 
of the gap lemma as stated in [Z.3], applied to the "lifted" matrix- valued ODE 



for the conjugating matrices P+. The x-derivative bounds < k < K -\- 1 
then follow from the ODE and its first K derivatives. Finally, the A-derivative 
bounds follow from standard interior estimates for analytic functions. □ 

Definition 2.2. Following lAGJf . we define the domain of consistent split- 
ting for the ODE system W = A{x, X)W as the (open) set of X such that the 
limiting matrix A_|_ is hyperbolic (has no center subspace ) and the boundary 
matrix B is full rank, with dim 5*+ = rankB. 



(2.5) 



W = (/ + 0+)Z 



P' = A+P - PA + (A - A+)P 
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Lemma 2.3. On any simply connected subset of the domain of consistent 
splitting, there exist analytic bases {f i, . . . , f fc}^ and {f fc+i, . . . , V]^}^ for the 
subspaces S+ and f/+ defined in Definition \2.2[ 

Proof. By spectral separation of f/+, S+, the associated (group) eigenprojec- 
tions are analytic. The existence of analytic bases then follows by a standard 
result of Kato; see jKat] . pp. 99-102. □ 

Corollary 2.4. By the Conjugation Lemma , on the domain of consistent 
splitting, the stable manifold of solutions decaying as x ^ +oo of (12.11) is 

(2.6) 5+:=span{P+t;+,...,P+i;+}, 

where W'\_ := P+Vj~ are analytic in A and C^^^ in x for A G . 

2.2 Definition of the Evans Function 

On any simply connected subset of the domain of consistent splitting, let 
Wi, . . . , Wf!' = P+Vi, . . . , P+v'^ be the analytic basis described in Corol- 
lary 12.41 of the subspace of solutions W of (12. ip satisfying the boundary 
condition W ^ aX +oo. Then, the Evans function for the ODE systems 
W = A{x, \)W associated with this choice of limiting bases is defined as the 
k X k Gramian determinant 



D{X) := det (mW^,...,MW^ 

(2.7) ; 

= det (MP+vt,...,MP+v 



|3;=0,A 

+ 
k 



|a;=0,A 



Remark 2.5. Note that D is independent of the choice of P+ as, by unique- 
ness of stable manifolds, the exterior products (minors) P+v^ A ■ ■ ■ A P+V^ 
are uniquely determined by their behavior as x ^ +oo. 

Proposition 2.6. Both the Evans function and the subspace are analytic 
on the entire simply connected subset of the domain of consistent splitting on 
which they are defined. Moreover, for A within this region, equation (12.11) 
admits a nontrivial solution W e L'^{x > 0) if and only if D{\) = 0. 



8 



Proof. Analyticity follows by uniqueness, and local analyticity of P+, f^. 
Noting that the first P+Vj^ are a basis for the stable manifold of fl2.1l) at 
X — > +CXO, we find that the determinant of BP+f^ vanishes if and only if B(A) 
has nontrivial kernel on 5+(A, 0), whence the second assertion follows. □ 

Remark 2.7. In this case that the ODE system describes an eigenvalue 
equation associated with an ordinary differential operator L, Proposition \2.6[ 
implies that eigenvalues of L agree in location with zeroes of D. (Indeed, they 
agree also in multiplicity; see fGJll \ GJ^ : Lemma 6.1, fZH^; or Proposition 
6.15 of IMZ^ .) 

When kerB has an analytic basis f . . . , f^.^, for example, in the 
commonly occurring case, as here, that B = constant, we have the following 
useful alternative formulation. This is the version that we will use in our 
analysis of the Green function and Resolvent kernel. 

Proposition 2.8. Let v^_^^., . . . .,v^^_f. he an analytic basis o/kerB, normal- 
ized so that det (B*, . . .v%) = 1. Then, the solutions Wj of (12. ip de- 
termined by initial data ^^^'(A, 0) = Vj are analytic in A and C^^^ in x, 
and 

(2.8) D{\) := det . . . , W.^' ■ ■ ■ ^W^) . 

Proof. Analyticity/smoothness follow by analytic/smooth dependence on ini- 
tial data/parameters. By the chosen normalization, and standard properties 

of Grammian determinants, -D(A) = det ( , . . . , , "^fc+i; ■ ■ ■ iV^n] i 
yielding ([23D- □ 

3 Construction of the Resolvent kernel 

In this section we construct the explicit form of the resolvent kernel, which 
is nothing more than the Green function G\{x, y) associated with the elliptic 
operator (L — A/), where 

(3.1) (L - A/)G'a(., y) = 6yl, Gx{0, y) = 0. 

Let A be the region of consistent splitting for L. It is an established fact 
(see [He]) that the resolvent (L — A/)~^ and the Green function Gx{x,y) are 
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meromorphic in A on A, with isolated poles of finite order. Gx in fact admits 
a meromorphic extension to a sector 

(3.2) ne = {X: Re{X)> -ei-e2\Im{X)\}, 61,62 >0. 

Writing the associated eigenvalue equation in the form of a first-order 
system (12.11) . we obtain 

(3.3) W' = A{X,x)W, MW{0) = 0, 
where 

'^'-t'j'^C- A=(^^_,;^,g_, AB-^-B-B-') 

and B = constant is the rank-n projection onto the first coordinate w of W, 
with kernel spanned by the constant basis v^_^_j = Cn+j, j = I, . . . ,n and Cj 
the jth standard basis element. 
Denote by 

(3.4) $o=(0O(^;A) ■•• 0O(x;A)) = (iy° ■■■ W^) 

and 
(3.5) 

<|.+ = (0+(x;A) •■■ 0+(x;A)) = (W-,+ 1 ■■■ W+) = {P+vt ■■■ P+v^) 

the matrices whose columns span the subspaces of solutions of (12. ip decaying 
at X = 0, +00 respectively, denoting (analytically chosen) complementary 
subspaces by 

(3.6) v^°=«(a;;A) ■■■ ^"(x; A)) = ■■■ W^) 
and 

(3.7) vl>+= (V;+(x;A) ■■■ V^+(x;A)) = ■■■ W+) . 

As described in the previous subsection, eigenfunctions decaying at both 

0, +00 occur precisely when the subspaces span$° and span$+ intersect, 

1. e., at zeros of the Evans function defined in (12.81) : 

(3.8) D^iX) := det($°, $+)|,.=o = (0? A ■ ■ ■ A 0° A 0+ A ■ ■ ■ A . 
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Lemma 3.1 ( |GZt IZHj ). For 61,62 > sufficiently small, Di is locally 
analytic on sector VLq as defined in (13.21) . 

Proof. Direct calculation showing that the domain A of consistent splitting 
is contained in VLq — B{0, r) for r > arbitrary and 6 sufficiently small, with 
extending analytically to -8(0, r). □ 

Lemma 3.2. Let H\{x,y) denote the Green function for the adjoint operator 
(L — A/)* on the half-plane x > 0. Then G\{x,y) = Hl{x,y). In particular, 
for X ^ y, the matrix z = Gx{x, .) satisfies 

(3.9) {z'B)' = -z'A + zX. 

Proof. Standard duality argument; see |ZHj for operators on the whole line. 

□ 

Considering (13.91) as an ODE system for the vector Z = {z, z'), it becomes 

(3.10) Z' = ZA{\,x), 
where 



(3.11) 



Afi-i - A'B-^ 



J -AB-^ - B'B-^ 

Lemma 3.3 ( |ZH] ). Z is a solution of (13.111) if and only if ZSW = constant 
for any solution W of (12. ip . where S = ^ . 

Proof. Direct computation/comparison with of {ZSW)'] see |ZHj . □ 
Using Lemma 13.31 we can define dual bases and by the relations 

(3.12) W^^^SW',^ = 6l 

Likewise, Ao,+ can be defined as 

We define also the dual subspaces 

(3.14) <l°=(0?(x;A) ■■■ 0°(x;A)) = (iy„\i ■■• Wl) , 
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(3.15) i.+ = (0+(a;;A) ■■■ 0+(x;A)) = ■•• , 

(3.16) ^°=(^?(x;A) ■■■ ^+(x;A)) = (lyo ■■■ W^) , 

(3.17) ^+=(V^+(x;A) ■■■ ^°(x;A)) = ■■■ W^2'n) • 

With these preparations, the construction of the Resolvent kernel goes 
exactly as in the construction performed in [ZHl IMZ3j on the whole line. 



Lemma 3.4. We have the the representation 
(3.18) 



Gx \ ^ I$+(A,x)M+(A)vl/0(A,y) for x > y, 

Gx^ Gx^J |<l>0(A,x)MO(A)^+(A,y) far x<y, 

where M°'+ are to be determined. 

Proof. See |ZH] Lemma 4.6. □ 

Using Lemma [3.41 we find the explicit coordinate-free representation for 
X > y: 

(3.19) ^'A = J'^-''U4z)S-\z)Ilo{z)T^-\ 
where 

(3.20) U+{y) = (<|.+ (y),0)($+(y),$-(y))"\ 

(3.22) J^^^" = (<l>+(x), $°(a:))(<l>+(z), $°(z))^\ 

(3.23) -p^-^y _ f n^A f ny) 

and similarly for x < y. 
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Corollary 3.5. The resolvent kernel may be expressed as 

.3 24^ G( . /(/n,0)<f+(x;A)M+(A)^o*(2/;A)(/„,0)*'- x > y, 
^' ^ ^ ^ \-(/„,0)$0(x;A)MO(A)^'+*(y;A)(/„,Or'- x<y, 

where 

(3.25) M(A) := diag(M+(A), M°(A)) = <l>~\z; X)S-\z)^-^* {z; A). 

4 Low-frequency bounds 

Our goal in this section is the estimation of the resolvent kernel in the critical 
regime |A| — >■ 0, i.e., the large time behavior of the Green function G, or global 
behavior in space and time. We are basically following the same treatment as 
that carried out for viscous shock waves of strictly parabolic conservation laws 
in |ZHt IMZ3j ; we refer to those references for details. In the low frequency 
case the behavior is essentially governed by the equation 

(4.1) Ut = L+U := -A+U^ + B+U^, 

Proposition 4.1. Assuming (H0)-(H4), let K he the order of the pole of Gx 
at X = and r be sufficiently small that there are no other poles in B{0, r). 
Then for X E such that |A| < r and for x > y > we have 

where djk{X) = 0{X~^) and ejk{X) = 0{X^~^) are scalar meromorphic func- 
tions, moreover K < order of vanishing of the Evans function D{X) at X = 0. 

Proof. See |ZH] Proposition 7.1 for the first statement and theorem 6.3 for 
the second statement linking order K of the pole to multiplicity of the zero 
of the Evans Function. □ 

Lemma 4.2. Assuming (H0)-(H4), for \X\ sufficiently small, the eigenvalue 
equation (L_|_ — X)W = associated with the limiting, constant- coefficient 
operator L+ has a basis of 2n solutions = e^^ ^"^^^VjlX) where fij' and 
Vj^ are analytic in X, consisting of n fast modes 

(4.3) ^;-7;+o(A). 

^ ' \// = 5+ + 0(A), 
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where 7^, are eigenvalues and associated right eigenvectors of B^-^A^, 
and n slow modes 

.4. /.^,,(A):=-A/a; + A^/5;/af + (9(A3), 

KX,(A):=r+ + 0(A), 

where , ^ , are eigenvalues and left and right eigenvectors of : = 
dF{u^) , and (3^ := l^B^r^^ > with i?+ := B{u^). The same is true for 
the adjoint eigenvalue equation 

(L+-A)*Z = 0, 

i. e, it has a basis of solutions 
with 

(4.5) v^+(A)=t;+ + o(A), 

(4.6) V:+,(A) = /; + 0(A), 
analytic in A. 

Proof See jMZ3] . □ 

Proposition 4.3. Assume (H0)-(H4) and ( 11. 7p . then, forr > sufficiently 
small, the Resolvent kernel G\ associated with the linearized evolution equa- 
tion 

(4.7) Ut = L+U := -A^U, + B+U,^ 

satisfies, for < y < x: 
(4.8) 

\d2d;Gx{x,t;y)\ < C(|Ar + e-'l"l)(|Ar + e"'!^!) ( ^ | e(-V4+AX+ | 



^{-X/a++\^l3+/af)x+iX/a+-\^P+/af)y\ 



+ E 

a+<0,a+>0 

< Iq^MtI < Ij ^ > 0, wzi/i similar bounds for < x < y. Moreover, each 
term in the summation on the righthand side of ( 14. Sp bounds a separately 
analytic function. 
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Proof. By 1.8 D does not vanish on Re{X) > 0, hence, by continuity, on 
|A| < r. Thus, according to (14.21) . all |(ijfc(A)| are uniformly bounded on 
|A| < r, and so it is enough to find estimates for fast and slow modes 0^, 0^, 
ipj' and ipj'. By using (13. 5p we find: 

and similarly for 0^, ^/^^ and ^z^^. Now using (12.41) and the fact, by (14. 4p . 

that e'^J^'^)^ is of order e"'"^' for fast modes and order e'^^"^ ^'^^^j +o(a^) 
for slow modes, substituting this and corresponding dual estimates in (14. 9 p 
and grouping terms, we obtain the result. □ 



5 High frequency bounds 

To analyze the high frequency behavior of the Green function of the boundary 
layer, we first establish some bounds for the projection terms in the Green 
function, using the symmetric formula 



(5.1) 



' G,{x,y) dyG,{x,y) \ ^ f ^^-^n+(x)5-^(y) tf x > y, 
^d^Gx{x,y) d^dyGx{x,y)J \:F'=^yil^{x)S-\x) if X < y. 



By setting x = \\2\x, A = j^, B{x) = w{x) = '"^(^) in the 

eigenvalue equation Lw = Xw associated with (II. 4p we obtain 

(5.2) W = MW + 0{\X-^\)W 

where 



(5.3) 



I 
\B 



and B' = (9(|A~2|) and |A| = 1. Since M{X,x) varies within a compact set, 
then there are eigenprojections Pq and P+ with property |P_^| = (9(|A^2 1) 
and |Pq| = (9(1 A~^ I) taking W onto the stable and unstable subspace. By 
using the two new coordinates = P^W and Yq = PqW, we obtain 

<-) a-)'^(-ir:)(^:)-<i-*i)e 
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Equivalently, we can find continuous invertible transformations Qo such 
that = Q^A^Q^^ and Eq = QoAqQq^ where 



(5.5) Re{E) := -{E+ + E^) < -p-^I. 



in the sense of quadratic forms. 

Again by coordinate change Z+ = Q+Y-^, Zq = QqYq we find 
where 

(5.7) ^ < kl < C\w\. 

From this we find by energy estimate that 

(5.8) (|z+|7<-2/3-i|^+|2 
and hence 

(5.9) < e-^^-^K^-^)! 

\z+{y)\ 

for any solution2;_|_ decaying at oo, where (3 < (3 and thus 

(5.10) \p^<e-^-h^-y^\ 

\z{y)\ 

for X > y, provided that |A| is sufficiently large. From this we obtain 

(,n) |M < c^e--'i<^-..i 

where C is as in ( 15. 7p . Applying a symmetric argument for the adjoint 
equation, we obtain the following lemma. 

Lemma 5.1. On the manifolds and defined in (13. 5p and (13.171) . for 

A sufficiently large, within the sector Qg = {\ : Re{\) > —6i — 92\Im{\)\} , 
61,62 > 0, we have in rescaled coordinates x, for some uniform C > 0, 



(5.12) IJ^^^^I, IJ^^^^I < Ce- 

for X > y and x < y respectively. 
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Lemma 5.2. In reseated coordinates x, A, for the projection terms n+(y) 
and n+(a;); the projection along $o onto for A sufficiently large, 

(5.13) |n+(y)|, \U4x)\<C 

for some uniform C > 0. 

fWi\ fW^ 
Proof. Choosing the coordinates I ^ 1 e C where — ' 



^^^1 I , we show 

Wo 



for small enough e and fixed c > such that 
projection along E — span {W^...W"-) onto F — span 



HI < ce and 



2, 

JH < c the 



(5.14) n := {w\..w'^){On, In){w\..w^^)-' 
satisfies 

(5.15) |n|<C2(c,e) 

To show this without loss of generality we assume that 

(5.16) {w-^\..w'-) = 1^^-^^ (w\..w-) = (^W^ 

Now it is sufficient to show that 

(5.17) \{w\...,w'-)\<C2{c,e). 
But, this amounts to showing that 



(5.18) 



M 4 

Jn O 

which amounts to showing that 



-1 



<C, 



(5.19) 



M /„ 
L O 



< C2{c), 



where \M\ < c. But, this is easy to show because 



and so 



M 7„ 
L. O 



M /„ 

In O 



In -M 

O L 



<l + \M\<C. 



□ 
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Proposition 5.3. Assume (H0)-(H4) and fll.7p . Then, for R > suffi- 
ciently large, the Resolvent kernel G\ associated with the linearized evolution 
equation (14.71) satisfies, for c, C > and 0<|a|,|7|<l.- 

(5.20) mG,{x^y)\ < C7|A|(^^)e-v^^ 

Proof. Recalling the coordinate-free representation (15.11) and combining with 
(I5.12P and (15.131) . we find that the Green function G\ in rescaled coordinates 
X, A satisfies 

(5.21) \d2d^G-^{xM<Ce-'-^, 

whence (I5.20p follows in the original coordinates. □ 

Remark 5.4. The argument of Lemma \5.2\ is the key new ingredient in the 
resolvent estimates for the boundary layer case as compared to the analysis on 
the whole line carried out for viscous shock layers in \ZH^ . making essential 
use of compatibility of the boundary condition with high-frequency behavior. 
On the whole line, there is no such requirement and high-frequency stability 
is automatic. 



6 Pointwise Green function bounds 

With the pointwise bounds established on the resolvent kernel G\, we obtain 
pointwise bounds on the Green function through the inverse Laplace trans- 
form formula by a simplified version of the stationary-phase arguments used 
in [ZH] for the shock case, repeated here for completeness. 

Proof of Theorem \1.3[ By sectoriality of L, we have the inverse Laplace trans- 
form representation (see [ZHj ): 

(6.1) G{t;x,y) = j^e^'Gx{x,y)dX. 

Let 6i > 0, 62 > be chosen sufficiently small, in particular so small as 
to satisfy the hypotheses of all previous assertions. By assumption (II. 7p . 
the large- 1 A| bounds on the Resolvent kernel, and analyticity of the Evans 
function Dl{X), it follows that Gx has finitely many poles in Qg (correspond- 
ing to roots of Dl), each with strictly negative real part. Choosing 61,62 
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still smaller, if necessary, we can thus arrange that G\ is analytic on ^Ig. It 
follows from Cauchy's Theorem that 



(6.2) G{x,t;y) = J^e^'Gx{x,y) dX, 

for any contour F that can be expressed as F = d{Qe \S) for S G C open. 

Case I. \x — y\/t large. We first treat the trivial case that \x — y\/t > S, 
S sufficiently large, the regime in which standard short-time parabolic theory 
applies. Set 

(6.3) «:=^Z^, R:=Pa', 

where (3 is as in (15. 5p . and consider again the representation of G, that is 

(6.4) G{x,t;y)= [ e^'Gx{x,y) dX, 

where Fi := dB{0, R)nCle and F2 := dne\B{0, R). Note that the intersection 
of F with the real axis is Xmm = R = PoP' . □ 

By the large |A| estimates of Proposition 5.3, we have for all A G Fi U F2 
that 

\Gx{x,y)\<C- 



Further, we have 
(6.5) 

ReX<R{l-T]u'^), XeTi,ReX < ReXo - r]{\ImX\ - \ImXo\), A G F2 

for R sufficiently large, where u is the argument of A and Aq and Aq are the 
two points of intersection of Fi and F2, for some rj > independent of d. 
Combining these estimates, we obtain 



(6.6) 



/ e^'GxdX\< [ C|A|-5e^^^*-^"^l^l"^l"-^lrfA 



+L 

<Ce-'^"'*/ R-h-'^^'^'^'' Rdu <Cr^e-^"'*. 
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Likewise, 

^Q^Re{\o)t-\i3\-^\M~^\^-y\ f |A|"^e(^'=^~-^'=^°)* MAI 

(6.7) Jr2 

<Ce-^^'* [ |/mA|-^e-^l-^"^^-^"*^°l* \dImX\ 

Jt2 

Combining these last two estimates, we have 

(6.8) \G{x,t;y)\ < Cr^e^e^^W^ < Cr^e-^^e^^W^, 

for 77 > independent of a. Observing that < < for any 

bounded a, for sufficiently large, we find that this contribution may be 

absorbed in any summand t~e mi . 

Case II. \x — y\/t bounded. We now turn to the critical case that \x — 
y\/t < S. A few remarks are in order at the outset. Our goal is to bound 
|G| by terms of form Cr^/^e-^'*/^, where a := (x - a+{t - \y/al\)/2t or 
a :— {x — y — alt)/2t are now uniformly hounded, by 

(6.9) |x-t/|/2t + max{|a+|}/2 < ^/2 + max |a+|/2. 

3 

Thus, in particular, contributions of order t~^^'^e~'^, rj > 0, can be ab- 
sorbed in any summand t~i/2e~(^~^~"fc /*^*) if wc take M sufficiently large. 
Likewise, for and Gy, contributions of order t~^e~^^ can be absorbed. We 
will use this observation repeatedly. 

In contrast to the previous case of large characteristic speed \x — y\/t > S, 
we are not trying to show rapid time-exponential decay. Rather, we are trying 
to show that the rate of exponential decay of the solution does not degrade 
too rapidly as a — 0: precisely, that it vanishes to order and no more. 
Thus, the crucial part of our analysis will be for small d. All other situations 
can be estimated crudely as described just above. 

Let r be sufficiently small that the small-| A| bounds hold on B(0, r) . Next, 
choose 61 and 62 still smaller than before, if necessary, so that Qg\B{0, r) C A. 
This implies that dflg fl B{0, r) ^ 0, giving the configuration pictured in the 
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Figure. Similarly as in the previous case, define F = Fi U where Fi is the 
portion of the circle dB{0,r) contained in fig, and F2 is the portion of dflg 
outside -8(0, r). 

(6.10) G{x,t;y)^ [ e'''Gx{x,y)dX+ [ e'''Gx{x,y) d\, 
we separately estimate the terms Jp^ and Jp^. 

Large- and medium-A estimates. The Jp^ term is straightforward. The 

points Ao, Ag where Fi meets F2 satisfy -Re(Ao) = —rj < 0. Moreover, com- 
bining the results low frequency case, we have the bound \Gx\ < C|A|~2 for 
A e F2. Thus, we have 

\[ e^^GxdX\<Ge-^^^°^ [ 1/^ Ar^e-"!-^"*^-^"*^"!* \dImX\ 

(6.11) Jr2 Jr2 

This contribution can be absorbed as described above. An analogous com- 
putation using IGa^tIj \G>^x\ — ^\M~^ shows that the F2 contribution to G^ 
and Gy is 0{t~^e~^*), and can likewise be absorbed. 

Small |A| estimates. It remains to estimate the critical term Jp^ e^^GxdX. 
This we will estimate in different ways, depending on the size of t. 

Bounded time. For t bounded, we can use the medium-A bounds \Gx\, 
\GxJ, \Gxy\ < C to obtain | Jp^ e^^GxdX\ < C2IF1I. This contribution is or- 
der Ce~'^ for bounded time, hence can be absorbed. 

Large time. For t large, wc must instead estimate Jp^ e^^GxdX using the 
small-|A| expansions. First, observe that, all coefficient functions djk{X) are 
uniformly bounded (since |A| is bounded in this case). 
Expanding G — J-^e^*Gx{x,y)dX as 

/ G GA ^ fxt f Gx Gx^ \ 

\Gy G:cyJ Jr \Gxy Gx^yJ 
we estimate the Jp^ contributions to G, Gx and Gy simultaneously. 
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Case II(i). {0 < y < x). By our low-frequency estimates, we have 



(6.12) 



where each rfjfc is analytic, hence bounded. We estimate separately each of 
the terms 

on the righthand side of fl6.12p . Estimates for terms 
go similarly. 

Case Il(ia). First, consider the critical case > 0, < . For this 

case, 

where 

' i^tiX) = -X/at + X'Pt/H f + 0{\') 
_p+(A) = -A/a; + \'Pp{aXf + ^(A^). 

Set 

a = , p := — > U. 

2t ' ^ t 

Define F'^^^ to be the portion contained in VLq of the hyperbola 
(6.13) 

Re{p^x-vty) + 0{\'){\x\ + \y\) 

= {l/at)Re[X{-atx/aj + y) + X\xf5j / {a^f - yp-J{atf)] 
= constant 

= {^/(h)[{Xmin{-alx/a^ + v) + Xl,i^{xl3^ al/ [a^f - yl3t/{alf )], 
where 

if |^|<e 



(6.14) Amin : = 
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Denoting by Ai, A^, the intersections of this hyperbola with dQg, define 
r'^ to be the union of AiAq and AqA^, and define T[ = U T[ . Note that 
X = a/p minimizes the left hand side of (16.131) for A real. Note also that that 
p is bounded for a sufficiently small, since a < e implies that 

{\atx/a+\ + \y\)/t<2\at\+2e 

i.e. + \y\)/t is controlled by a. 

With these definitions, we readily obtain that 

ReiXt + p+x - uty) < -(t/a, )(aV4p) - r]Im{X)H 
< -ah/M -r]Im{X)H, 

for A G r'^^ (note: here, we have used the crucial fact that a controls (|x| + 
\y\)/t, in bounding the error term (9(A^)(|x| + \y\)/t arising from expansion 
Likewise, we obtain for any q that 

(6.16) / lAI-^e^^^^^+^^-'^-^^c/A < Cr^-ie-"'*/^^ 



for suitably large C, M > (depending on q). Observing that 

« = (4/4)(^-4(^-l2//4l))/2t, 

we find that the contribution of (I6.16P can be absorbed in the described 
bounds for t > \y/a^\. At the same time, we find that a > x > for 
t < Izz/'^fc I) whence 

a>{x-y- ap)/Mt + \x\/M, 

for some e > sufficiently small and M > sufficiently large. 
This gives 

provided \x\lt > a'^, a contribution which can again be absorbed. On the 
other hand, if t < |x/a^|, we can use the dual estimate 

.g^^. (^ = i-y-4it~\x/a+\))/2t 

> {x-y-att)/Mt+\y\/M, 
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together with \y\ > to obtain 

g-"Vp < g-(2^-3/-4*)V^-fte-'?|y|^ 

a contribution that can hkewise be absorbed. 

Case Il(ib). In case < or > 0, terms \(pj'\ < Ce~^^^^ and 
\ip^\ < Ce"'''^' are strictly smaller than those already treated in Case Il(ia), 
so may be absorbed in previous terms. 

Case Il(ii) {0 < x < y). The case < x < y can be treated very similarly 
to the previous one; see |ZH] for details. This completes the proof of Case 
II, and the theorem. 



7 Nonlinear Analysis 

Introducing the perturbation variable 

(7.1) u{x,t) := u{x,t) — U{x), 
we obtain 

(7.2) ut- Lu = Q{u)^, 
where the second-order Taylor remainder satisfies 

(7.3) Q{u) := f{u + u)~ f{u) - df{u)u = 0{\u\^) 
so long as |m| remains bounded. 

Lemma 7.1 (Integral formulation). Under the assumptions of Theorem \1.5[ 
there exists a classical solution of fl7.2p for 0<t<T,T>0, continuous in 
L°°{x) at t = 0, extending for all t > such that u{-, t) remains sufficiently 
small in fl given by 



u{x,t) = I G{x,t;y)g{y)dy 

(7.4) 



Jo 



[ Gy{x,t - s;0)Bh{s)ds 
Jo 

Gy{x,t - s; y)Q{u){y, s) dyds. 
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Proof. From Lemma fl3.2l) and the inverse Laplace representation fl6.1l) we 
find that G{x,t — s;y) considered as a function of y,s satisfies the adjoint 
equation 

(7.5) ids-LyyG*{x,t--,-) = 0, 
or 

(7.6) - - {GA)y + GAy = {GyB)y. 

Likewise, reviewing the construction of the resolvent, we find G\{x,Q) = 0, 
yielding 

(7.7) G(x,t- s;0) = 0. 

That is, G*{x,t — ■; ■) is the Green function for the adjoint equation, as may 
alternatively be seen directly by a duality argument analogous to the proof 
of Lemma (13.21) . 

Thus, integrating G against (17. 2p . integrating by parts, and using the fact 
that G = and u = h on the boundary y = 0, we obtain for any classical 
solution of (17. 2p that 
(7.8) 



JO 



G{x,t- s; y)Q{u{y, s))ydyds = 

G{x, t — s;y){ds — Ly)u{y, s) dy ds 

/ / {{ds- LyyG*y{x,t- s;y)u{y,s)dyds 
Jo Jo 



Jo 



+ u{x,t)- G{x,t;y)g{y)dy 



[ Gy{x,t - s;0)Bh{s)ds, 
Jo 



from which we obtain (17. 4p by rearranging and integrating by parts the term 

Indeed, (17.40 may be taken as the definition of a weak solution in L°°{x, t). 
(One can see using convolution identities that this agrees with the usual 
definition in terms of integration against test functions G C^(]R x M).) 
Existence of weak solutions can be obtained by a standard contraction map- 
ping/continuation argument using the convolution bounds of Lemmas 17. 2117^ 
below; we omit the details, since we shall carry out quite similar but more 
difficult estimates in the proof of stability. Smoothness of solutions may then 
be obtained by a bootstrapping argument as sketched in Appendix \^ □ 
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To establish stability, we use the following lemmas proved in [HZj . 
Lemma 7.2 (Linear estimates |HZ] ). Under the assumptions of Theorem 

m 

r+oo 

(7.9) / \G{x, t- y)\{l + dy < C{e + + ^2)(x, t), 



for < t < +00, some C > 0. 

Lemma 7.3 (Nonlinear estimates [HZj ). Under the assumptions of Theorem 

ft f + CO 

(7.10) / / \Gy{x,t-s-y)\m{y,s)dyds<C{e + ^i+ij2){x,t), 



JO JO 

for <t < +00, some C > 0, where 

,^ *(l/, s) := (1 + sf/'s-"\e + ^1 + i,,f{y, s) 

+ {i + sr\e + ij^ + ,lj2){y,s). 

We require also the following estimate accounting boundary effects. 
Lemma 7.4 (Boundary estimate). Under the assumptions of Theorem \1.5[ 

(7.12) f Gy{x,t - s;0)Bh{s)ds < CEq{6 + + ^2){x,t), 

JO 

for < t < +00, some C > 0. 

Proof. The estimate on ^, where Gy{x,t — s; 0) is nonsingular, follows 
readily by estimates similar to but somewhat simpler than those of Lemma 
07.31) . which we therefore omit. 

To bound the singular part J^_-^^, we integrate (17. 6p in y from to +oo, 
recalling that G{x, t — s;0) = 0, to obtain 

P + OO P + OO 

(7.13) GyB = - Ay{y)Gix, t-s;y)dy- Gs{x, t-s;y) dy. 

JO JO 



26 



Substituting in the lefthand side of fl7.12p . and integrating by parts in s, we 
obtain 



GyBh{s)ds = i^j Ay{y)G{x,T]y)dy^h{t-T)dT 

(7.14) - ( G{x, T- y) dyY'{t - r) dr 

"+00 



+ G{x,l;y)dyy{t-1), 



which by J \G\dy < G has norm bounded by maxo<r<i(|/i| + \h'\){t — r) 
Combining this with the more straightforward estimate 



t-i 



Gy{x,t;0)Bh{s)ds < [ \Gy{x,T;0)\Bh{s) ds 

Jo 

< G max \h{t - t)\ C t-^-^''^"'^^ dr 

0<r<l Jq 



(7.15) =C|x|-2 max |/i(t-r)| 

0<T<1 ' 

X l\\x\^/T)e-\^\'"'^dT 



< G max \h{t — r)||a;| 

0<T<1 



we find that the contribution from J^*_^ has norm bounded by 



max(|/i| + |/i'|)(t-r)(l + |x| 

0<r<l 



\-2 



Combining this estimate with the one for /g , we obtain (17J2|) . □ 

With these preparations, the proof of stabihty is straightforward. 
Proof of Theorem Define 

(7.16) at)-= sup |n|(^ + ^i + ^2)"'(2/,t). 

y,Q<s<t 

We will establish: 
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Claim. For all t > for which a solution exists with ( uniformly bounded 
by some fixed, sufficiently small constant, there holds 

(7.17) at)<c2{Eo+aty). 

From this result, provided Eq < 1/4C|, we have that ({t) < 2C2-E0 implies 
((t) < 2C2-E05 and so we may conclude by continuous induction that 

(7.18) C{t)<2C2Eo 

for all t>0. (By Lemma [7.11 and standard short-time estimates, u G C^{x) 
exists and ( remains continuous so long as ( remains bounded by some uni- 
form constant, hence fl7.18p is an open condition. From fl7.18p and the defi- 
nition of C in fl7.16p we then obtain the bounds of (11.141) . Thus, it remains 
only to establish the claim above. 

Proof of Claim. We must show that u{9 + ipi + 4'2)~^ is bounded by 
C{EQ + ((t)^), for some C > 0, all < s < t, so long as ( remains sufficiently 
small. By fl7.16p . we have for alH > and some C > that 

(7.19) \u{x,t)\<Cit)i9 + iP, + ^2){x,t), 
and therefore 

(7.20) \Q{u){y,s)\<Caty^{y,s) 

with \E' as defined in (17.111) . for < s < t. Combining (I7.20p with represen- 
tation (17.41) and applying Lemmas l7.2H7.4t we obtain 

\uix,t)\< \Gix,t;y)\\g{y)\dy+ Gy{x,t - s;0)Bh{s) ds 











(7.21) 



pt poo 

+ \Gy{x,t - s;y)\\{Q{u)){y,s)\dyds 

Jo Jo 

00 

-3/2, 



<Eo / \Gix,t;y)\il + \y\)-'^'dy 
Jo 

+ / Gy{x,t- s;0)Bh{s)ds 
Jo 

pt poo 

+ GC{tf / \Gy{x,t-s-y)\m{y,s)dyds 
Jo Jo 

<G{E, + atf){e + i,^+^2){x,t). 

Dividing by (6'+'?/'i+'02)(a;, ^), we obtain (I7.17P as claimed. This completes 
the proof of the claim, and the theorem. □ 
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A Smoothness of solutions 

In this appendix, we briefly sketch the proof that weak solutions defined by 
(17.41) are necessarily smooth, classical solutions as well, by indicating how to 
get the necessary derivative bounds. 

Time-derivative. Rewriting the second, boundary term, on the right- 
hand side of (17.41) using its convolution structure, as 

t 

Gy{x,T] 0)Bh{t-T) dr, 

and differentiating in t, we obtain 

Gy{x,t;0)Bh{0) + [ Gy{x,T;0)Bh'{t-T)dT, 
Jo 

for which the first term is bounded and smooth for x,t > 0, and the second 
by the same estimate as in (I7.15p is bounded by 



C\x\ 



[ \h'{t-T)\dT <G\x\-Hog{l+t). 

Jo 



Differentiating the first and third terms with respect to t and integrating 
the third term by parts in y yields 



(A.l) 



POO 

G{x,t;y)g{y)dy- I / Gy{x,t - s;y)Q{u)s{y, s) dyds 



t/2 Jo 



t/2 foo 

Gyt{x, t- s; y)Q{u){y, s) dyds, 



Jo 



from which, in combination with the boundary estimate already performed, 
we may readily obtain a short-time bound \ut\ < C\x\^'^t^^ by Picard itera- 
tion. 

Spatial-derivatives. Likewise, differentiating (I7.14p with respect to x, 
we may bound the x-derivative of the boundary term GyB ds by 

C fr-'^Wh'l + \h\){t -r)\dT< Clog(l + t). 
Jo 
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Differentiating the first and third terms of the righthand side of fl7.4p with 
respect to x and integrating the third term by parts in y yields 



ft /"OO 

G{x,t;y)g{y)dy- I / G^{x,t ~ s]y)Q{u)y{y, s) dyds 



(A.2) 



t/2 Jo 



t/2 poo 

Gy^{x,t- s] y)Q{u){y, s) dyds, 



JO 



from which, in combination with the boundary estimate already performed, 
we obtain a short-time bound \ux\ < Ct~^^'^ by Picard iteration. From the 
bounds on \ut\ and \ux\, finally, we obtain bounds on lu^^l by the equation 
satisfied by u. 
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